The spherical 2 + p spin glass model: an exactly solvable model for glass to spin-glass 

transition. 
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We present the full phase diagram of the spherical 2 + p spin glass model with p > 4. The main 
outcome is the presence of a new phase with both properties of Full Replica Symmetry Break- 
ing (FRSB) phases of discrete models, e.g, the Sherrington-Kirkpatrick model, and those of One 
Replica Symmetry Breaking (IRSB). The phase, which separates a IRSB phase from FRSB phase, 
is described by an order parameter function q{x) with a continuous part (FRSB) for x < m and a 
discontinuous jump (IRSB) at a; = m. This phase has a finite complexity which leads to different 
dynamic and static properties. 
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In the last years many efforts have been devoted to the 
understanding of complex systems such as spin glasses, 
structural glasses and others. Common denominator 
of all these systems is a large number of stable and 
metastable states whose complex structure determines 
the their static or dynamic behaviors. In this framework 
mean-field models, and among them spherical models, 
represent a valuable tool of analytical and theoretical 
investigation since they can be largely solved. Up to 
now only spherical models with One Replica Symme- 
try Breaking (IRSB) phases were studied, mainly due 
to their relevance for the fragile glass transition 0,M Q • 

To our knowledge the possibility infinite or Full Replica 
Symmetry Breaking (FRSB) phases in spherical models 
was first pointed out by Nieuwenhuizen on the ba- 
sis of the similarity between the replica free energy of 
some spherical models with multi-spin interactions and 
the relevant part of the free energy of the Sherrington- 
Kirkpatrick (SK) model HQ. In this paper Nieuwen- 
huizen presented some results for the FRSB phase but 
a complete analysis was not provided. The problem was 
considered some years later in connection with the 
possible different scenarios for the critical dynamics near 
the glass transition j^. This paper, however, analyzed 
only the dynamical behavior in the IRSB phase. 

The Model. The model we shall consider is the spher- 
ical 2+p spin glass model without external field defined 
by the Hamiltonian 
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where J- / , are uncorrelated zero mean Gaussian vari- 
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and Ui are N continuous variables obeying the spheri- 
cal constraint ^ ■ of — N . The properties of the model 



strongly depend on the value of p: for p — i the model 
reduces to the usual spherical p-spin spin glass model 
in a field Q with a low temperature IRSB phase, while 
for p > 4 the model posses an additional FRSB low- 
temperature phase [3- A partial analysis of the phase 
space of the model was carried out in Ref. where, how- 
ever, only the dynamical stability of the IRSB phase was 
considered leaving out a large part of the phase space 
and in particular the question of the transition between 
the IRSB and the FRSB phases. 

In this Letter we complete the study of the phase space 
focusing in particular on the transitions lines. We have 
studied the model by three complementary approaches. 
The first employs the replica method and analyze the 
disorder-averaged logarithm of the partition function fol- 
lowing Ref. U The second approach stars form the mi- 
croscopic dynamics and extend the results of Ref. ^ while 
the latter uses the Thouless- Anderson-Palmer approach 
[loj . In this Letter we shall mainly follow the replica 
approach, discussing differences with other approaches 
when necessary. A complete analysis of the properties of 
the model is beyond the scope of the Letter and will be 
presented elsewhere. 

Applying the standard replica method the free energy 
per spin / can be written as function of the symmetric 
n X n replica overlap matrix Qafs as 



-pf = -Pfo + s{oo) + lim - max G[Q] (3) 

On Q 

where /o is an irrelevant constant, s(oo) = (14- ln27r)/2 
the entropy per spin at infinite temperature T ~ 1//3, 
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where we have used the shorthand fip = {pj'^^^Y /2p. 
The spherical constraint is ensured by the condition 

Q aa Q 1 ■ 

Following Parisi 11] the overlap matrix QafS for a num- 
ber R of steps in the replica symmetry breaking is divided 
into successive boxes of decreasing size Pr, with pq — n 
and pr+1 = 1. The replica symmetric case and the FRSB 
case are obtained for i? = and i? — *■ oo, respectively. In 
the Parisi scheme the elements of Qaf3 are then given by 
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with Qr+i = q, where the notation aOP = r means that 
a and (3 belong to the same box of size Pr but to two 
distinct boxes of size Pr+i < Pr- The matrix obtained is 
conveniently expressed using the function 
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which equals the fraction of pair of replicas with overlap 
less or equal to q. Inserting this structure into eqs. ©- 
neglecting terms of order O(n^), and replacing the 
sums by integrals, one gets after a little of algebra, 



213 f - 2s(oo) - 2/3/0 + / dqx{q) —g{q) 
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where q{l) = qn and q{x) is the inverse of x{q). Max- 
imization of / with respect to q{x) leads to the self- 
consistent equation(s) for the order parameter function 
q{x). Depending on the value of the coupling strengths 
J^P^ and of the temperature T the function q{x) displays 
different forms which characterize the different phases of 
the model. Figurenshows the phase diagram in the space 
of the "natural" parameters — fj,2 for p = A. In the 
following we shall limit ourself to the case p = 4, how- 
ever the results are qualitatively valid for any p > 4. The 
analysis of the figure reveals four different phases, which 
will be discussed in the forthcoming part of this Letter. 

The Paramagnetic phase (PM). This phase exists for 
not to large values of coupling parameters strengths 
and/or high temperature and is characterized by a null 
order parameter function: q{x) = in the whole range 
X e [0, 1]. The phases becomes unstable above the line 
/i2 = 1 (DeAlmeida-Thouless line) where the "replicon" 
A = 1 — /i2 becomes negative. In this region for p > 4 and 
fip not too large a IRSB solution is also unstable and a 
more complex phase (FRSB) appears. Below fi2 — ^ the 
PM phase remains stable for all values of ^ip, similarly 
to what happens in the spherical p-spin model without a 
field 0], however as /Xp increases a more thermodynam- 
ically favorable IRSB phase with a non vanishing order 
parameter appears. 
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FIG. 1: Phase digram of the spherical 2 + 4 spin glass 
model. PM; Paramagnetic phase; IRSB: one-replica sym- 
metry breaking phase; FRSB: full replica symmetry break- 
ing phase; l-FRSB: one-full replica symmetry breaking. The 
dashed lines refers to dynamics. The m-line shown are ob- 
tained for m = 0.5. The continuous transition between the 
PM and the FRSB phases and between the FRSB and 1- 
FRSB phases are the same for statics and dynamics. Inset: 
the discontinuous transition between the FRSB and the 1- 
FRSB phases. 



The One Replica Symmetry Breaking phase (IRSB). 
This phase is characterized b y a step- like order parameter 
function q{x) = qi9{x — m) |l2| and is stable as long as 
the replicon eigenvalue is positive: 
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Maximization of / with respect to qi and m leads to 
the IRSB equations whose solution can be conveniently 
expressed defining qi — (1 — t/)/(l — y + my) in term of 
the function 



z{y) = -2y 
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introduced in Ref. 3 for the solution of the spherical 
p-spin spin glass model. For p = 4 the solution reads 
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By fixing the value of m G [0, 1] and varying y from j/min : 

z(ymin) = 1/2 (a<2 = 0) tO?/,nax : ■^(ymax) = l/2(l-|-2/i„ax), 

where the replicon Q vanishes, one obtains the so called 
TO-lines. The transition between the PM and the IRSB 
phases corresponds to the m — 1-line. Along this line qi 
jumps discontinuously from zero (PM) to a finite value 
(IRSB) however, since m = 1, the thermodynamic quan- 
tities remain continuous. Inserting into (|ll|l the value 
ymax, for which the replicon vanishes, and varying m from 
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FIG. 2: Schematic form of the order parameter function q{x) 
in the 1-FRSB phase. 



1 to one obtains the critical line between the IRSB and 
the 1-FRSB phase. 

The static approach requires that / be maximal with 
respect to variations of m. The dynamics, on the other 
hand, leads to the different conditions (marginal condi- 
tion) 



(12) 



which can be stated by saying that the derivative off 
with respect to m (the complexity) be maximal |l3l Il4l 
lisl . As a consequence the transition lines for dynam- 
ics and statics do not coincide. Due to space limitations 
the equations for the dynamical transition lines will not 
be reported, but only drawn in Figure^for completeness. 

The One-Full Replica Symmetry Breaking phase (1- 
FRSB). The analysis of the instabihty of the IRSB solu- 
tion reveals that in order to stabilize the phase above the 
line where the replicon Q vanishes a non-zero qo would 
be needed. However in the absence of external fields the 
order parameter function must vanish as a; — > 0, and 
hence a IRSB solution is not possible. On the other 
the different location of the static and dynamic instabil- 
ity lines suggests that sonic sort of IRSB-like form must 
survive in the solution. The way out is to look for a so- 
lution which below go has a structure which vanishes as 
a; — > 0. The most general form is an order parameter 
that has a discontinuity at x — m, is continuous below it 
and vanish for a; = 0: 
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for X > m 



qo{x) for X < m 



(13) 



with q{0) = and hm^_„- q{x) = (?o 7^ liTa^^jn+ qix) = 
q{l) = qi, see Figure |21 That this is the correct ansatz 
also follows from the numerical solution of the Parisi 
equations derived from stationarity of / with respect to 



The 1-FRSB equations are obtained by inserting the 
form H13() into the replica free energy (jSJ and imposing 
stationarity with respect to qo{x), qi and m. The result- 
ing equations can be solved in term of m-lines similarly 
to what done for the IRSB case. For the p = 4 case the 
solution for the "discontinuous" part of q{x) reads 
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and 51 = (1 - y)/[l - y 
part of q{x) is given by 



my(l — t)]. The "continuous" 



q = dq'[fi2 + 3^ 
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Xiq) = l-qi+ m{qi - qo) + / dq'x{q') (17) 

Jq 

For any fixed value of m G [0,1] these equations can be 
solved varying y from y,nin ■ x{y,nin) — (transition line 
to IRSB phase) up to y = 1 {t{y = 1) = 1) where the 
difference between qi and qi vanishes. These lines are 
the continuation into the 1-FRSB phase of the IRSB m- 
lines. In particular the m = 1 line represents the transi- 
tion between the 1-FRSB to the FRSB phase. For many 
aspects this transition is similar to the transition between 
the PM and the IRSB phases, indeed qi jumps discon- 
tinuously from a null value (FRSB) to a finite value (1- 
FRSB) however the discontinuity appears at m = 1 so 
that the thermodynamic quantities are continuous across 
the transition. The critical m = 1-line ends at the point 
where qi — qo, which for p = 4 is 



qi = qo 
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From this end-point on the transition between the FRSB 
and the 1-FRSB phases can only take place without a 
jump in in order parameter function. The value of t, the 
ratio between qi and go, increases along the m-lines as 
one moves away from the transition line with the IRSB 
phase, and the lines terminates when t = 1 {qo = qi). 
The set of all end-points for to S [0,1] defines the con- 
tinuous critical line between the 1-FRSB and the FRSB 
phases, which for p = 4 reads 
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On this line qi = qo = 1/(1 + 3m) and Xc = x{qQ) = 
m. It is worth to note that in going from the 1-FRSB 
to the FRSB phase the solution changes from stable to 
marginally stable, and remains marginally stable in the 
whole FRSB phase. 

The presence of a discontinuity in in the order parame- 
ter function leads to a finite complexity so that static and 
dynamics calculations leads to different solutions being 
the first associated with states of smallest (zero) com- 
plexity while the latter with states of largest complexity. 
As a consequence the m-lines in the two cases are differ- 
ent, as shown in Figure^ We shall not report the expres- 
sion for the dynamical m-line, this will be done elsewhere. 
The inset of the figure shows the different transition lines 
between the FRSB and the 1-FRSB phases. The discon- 
tinuity, and hence the complexity, vanishes on the con- 
tinuous transition line and two solutions coincide on this 
line and in the whole FRSB phase. 

The Full Replica Symmetry Breaking phase (FRSB). 
In this phase the order parameter function is continuous 
and given by eq. H16|l with qi = qo- By expressing (|16|l in 
term of q{x), instead of x{q), ad taking successive deriva- 
tives whit respect to x a power expansion of q{x) about 
X — can be computed |17l Il9l | . For the 2-1-4 model it 
turns out that q{x) contains only odd powers of x, the 
first of which are 

3/2 7/2 ,„ 11/2 

q{x) ^l^a: + l^x^ + i^gV + (20) 

Using this expression one can show that as the PM-FRSB 
transition line is approached from above r = /i2 — 1 ^ 0^ 
then both qo = q{xc) and Xc vanishes linearly with r as 

<1o = 1+0{t% a;, = ^T + 0(r2), r ^ 0+. (21) 

so that the transition between the PM and the FRSB 
phases takes place continuously without any jump in the 
order parameter function. 

Conclusions. To summarize in this Letter we have pro- 
vided the full phase diagram of the spherical 2 + p spin 
glass model with p > 4. Despite its simplicity the model 
has a rather rich phase. Indeed not only it presents a 
IRSB phase similar to that of the spherical p-spin spin 
glass model and a FRSB phase similar to that of the SK 
model, but it also posses a completely new phase with an 
order parameter made of a continuous part for x < m < I 
much alike the FRSB order parameter and a discontin- 
uous jump at X — m resembling the IRSB case. To 
emphasize its mixed nature this phase, which separates 
the FRSB phase and the IRSB phase, has been called 
1-FRSB. For may aspects the 1-FRSB phase is similar 
to the IRSB, in particular it is stable, at difference with 
discrete models as the SK model where a solution of this 
form is always unstable and has a finite complexity 
counting metastable states that are strict minima of the 
free energy landscape. 



The transition between the FRSB and the 1-FRSB 
phase can be either continuous (for /Z2 large enough) or 
discontinuous. In the first case the transition line is the 
continuation of the discontinuous transition between the 
PM and the IRSB phases and as for the latter the discon- 
tinuity appears at m = 1 . The presence of finite complex- 
ity in the 1-FRSB phase makes the static and dynamics 
transition different. The two lines join together at the 
end-point where the discontinuity at m = 1 in the or- 
der parameter function vanishes. From this point on the 
transition between the FRSB phase and the 1-FRSB can 
only take place continuously with the discontinuity of the 
1-FRSB phase which vanishes at the transition. Along 
the the continuous transition line the complexity vanishes 
and the static and dynamic approaches lead to the same 
results, in agreement with the conjecture made in Ref. 
i21i that the FRSB phase has vanishing complexity. 

In conclusion we believe that this is a rather promising 
model since not only it can be fully solved, but it possess 
different phases which can be fully analyzed. Moreover 
it posses an interesting transition between two different 
glassy phases, similar to what found in some colloidal 
suspensions |2a |. 
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